function [c,ceq] = reflection_con(v)
% Nonlinear inequality constraints
c = [];
% Nonlinear equality constraints
ceq = v(8)*v(5)/v(10)-v(3);
end